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Abstract. We present the galaxy rest-frame near-IR Luminosity Function (LF) and its cosmic evolution to 2; ~ 1.5 
based on a spectroscopic survey of a magnitude limited sample of galaxies with Ks < 20 (the K20 survey, Cimatti 
et al. 2002b) . The LFs have been derived in the rest-frame J and Ks bands. Their evolution is traced using three 
different redshift bins {zmean — 0.5, 1, 1.5) and comparing them to the Local near-IR Luminosity Function. The 
luminosity functions at different redshifts are fairly well fitted by Schechter functions at 2; < 1.3. The faint-end 
of the LFs {L < L*) is consistent with the local estimates, with no evidence for a change either in the slope or 
normalization up to z < 1.3. At higher redshift this part of the luminosity function is not well sampled by our 
data. Viceversa, the density of luminous galaxies {Mk^ ~ 5 log/iro < —25.5) is higher than locally at all redshifts 
and relatively constant or mildly increasing with redshift within our sample. The data are consistent with a mild 
luminosity evolution both in the J- and -f^s-band up to z ~ 1.5, with an amplitude of about AA/j ~ —0.69 ±0.12 
and AM/f ~ —0.54 ± 0.12 at z ~ 1. Pure density evolution is not consistent with the observed LF at z < 1. 
Moreover, we find that red and early-type galaxies dominate the bright-end of the LF, and that their number 
density shows at most a small decrease (< 30%) up to z ~ 1, thus suggesting that massive elliptical galaxies were 
already in place at z ~ 1 and they should have formed their stars and assembled their mass at higher redshift. There 
appears to be a correlation of the optical/near-IR colors with near-IR luminosities, the most luminous/massive 
galaxies being red/old, the low- luminous galaxies being instead dominated by blue young stellar populations. We 
also investigate the evolution of the near-IR comoving luminosity density to z ~ 1.5, finding a slow evolution with 
redshift (px{z) = p\{z = 0)(l + z)'''^^ with /?(J) ~ 0.70 and l3{Ks) ~ 0.37). Finally, we compare the observed LFs 
with the predictions of a set of the most updated hierarchical merging models. Such a comparison shows that the 
current versions of hierarchical models overpredict significantly the density of low luminosity galaxies at 2 < 1 
and underpredict the density of luminous galaxies at 2 > 1, whereas passive evolution models are more consistent 
with the data up to 2 ~ 1.5. The GIF model (Kaufmann et al. 1999) shows a clear deficiency of red luminous 
galaxies at z ^ 1 compared to our observations and predicts a decrease of luminous galaxies with redshift not 
observed in our sample. 
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- infrared: galaxies 
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1. Introduction 

Over the past few years, a wealth of observations from 
deep surveys of optically-selected high-redshift galaxies 
(e.g. Madau et al. 1996, Steidel et al. 1999), complemented 
by observations in the far-IR/sub-mm (Hughes et al. 1998, 
Barger et al. 1999), allowed significant progress in our un- 
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derstanding of the evolution of galaxies from the present- 
epoch back to z ~ 4 and beyond (Steidel et al 1999; 
Madau, Pozzetti, & Dickinson 1998). However, since these 
samples of high- 2; galaxies were all selected at optical or 
far-IR/sub-mm wavelengths, they are dominated by ob- 
jects with on-going star formation. Therefore, such stud- 
ies placed constraints more on the evolution of the star 
birth rate activity than on the formation and assembly of 
stellar systems through cosmic time. 

The study of faint galaxy samples selected in the near- 
infrared represents an important and complementary pos- 
sibility to address the still open questions on how the 
formation and evolution of massive systems evolved with 
time compared to the predictions of the different theoreti- 
cal scenarios (Broadhurst et al. 1992). Another advantage 
of the near-IR selection (in particular in the _fC-band) is 
that the k-corrections are relatively insensitive to galaxy 
type and fairly small also at high redshift (Cowie et al. 
1994), and the dust extinction effects are less severe than 
in optical samples. 

Since the rest-frame near-IR light is a relatively good 
tracer of the galaxy stellar mass (Gavazzi et al. 1996; 
Madau, Pozzetti & Dickinson 1998) the near-IR galaxy 
Luminosity Function (LF) can provide a reasonable esti- 
mate of the Galaxy Stellar Mass Function (GSMF). Only 
recently the 2MASS (Jarrett et al. 2000) surveys allowed 
accurate determinations of the local near-IR luminosity 
and of the Galaxy Stellar-mass functions (Cole et al. 2001, 
Kochanek et al. 2001), while only few attempts have been 
made to reconstruct their evolution with redshift using 
deep surveys of near-IR selected samples (i.e. Cowie et al. 
1996; Cohen et al. 1999, Cohen 2002). 

In order to address the above questions, we per- 
formed a new spectroscopic survey of a complete sam- 
ple of galaxies selected with Kg < 20 (the K20 sur- 
vey; http : //www . arcetri . astro . it/~ k20/). The sur- 
vey and the sample are described in detail in Cimatti et al. 
(2002b, hereafter Paper III), while the spectral and clus- 
tering properties of the Extremely Red Objects (EROs) 
are discussed in Cimatti et al. (2002a, Paper I) and Daddi 
et al. (2002, Paper II), respectively. The redshift distri- 
bution for the whole sample of Kg-hand selected galaxies 
is given in Cimatti et al. (2002c, Paper IV). Here we re- 
call that the K20 sample includes 546 objects to Kg < 20 
(Vega system), selected from a 32.2 arcmin^ area of the 
Chandra Deep Field South (CDFS; Giacconi et al. 2001) 
and from a 19.8 arcmin^ field centered at 0055-269. The 
total area of the two fields is 52 arcmin^. Optical multi- 
object spectroscopy was mainly obtained with the ESO 
VLT -I- FORSl and F0RS2, while a small fraction (~ 4%) 
of the objects was observed with near-IR spectroscopy us- 
ing the VLT + ISAAC. We have imaged the 0055-269 
field over 10 bands {U BGV RRwIzJKg), obtained with 
the ESO NTT + SUSI2 (UBGVR^I) and SOFI {J Kg), 
and VLT FORSl (i? and z), for a total of about 45 
hours of integration. In the CDFS field, we used a combi- 
nation of pubhc EIS NTT data {UKg) and deep FORSl 
images {BVRIz, courtesy of P. Rosati & M. Nonino). 



The spectroscopic redshift completeness is 94% and 87% 
for Kg < 19 and Kg < 20, respectively, and it increases 
to 98% if we include the photometric redshifts obtained 
with the deep multi-band imaging for the spectroscopi- 
cally unidentified or unobserved objects (Cimatti et al. 
2002b). The K20 sample is the largest and most complete 
spectroscopic sample of galaxies with Kg < 20 available 
to date. 




redshift 

Fig. 1. The magnitude- redshift diagram for all the galaxies in 
the spectroscopic sample to Ks < 20. Filled circles represent 
galaxies identified spectroscopically, while empty squares are 
unidentified or unobserved galaxies plotted at z = Zphot- 



In this paper, we investigate the evolution of the near- 
IR luminosity function up to z ~ 1.5 based on the K20 sur- 
vey sample. Thanks to the higher statistical significance 
and completeness of our sample, it is possible to use the LF 
evolution to place new and more stringent constraints on 
the formation and evolution of massive galaxies than was 
possible from previous surveys (Cowie et al. 1996, Cohen 
et al. 1999, Cohen 2002). In particular our survey has the 
advantage to include a complete sample of EROs, partially 
with spectroscopic identifications, which usually were not 
included in previous spectroscopic surveys (e.g. Cohen 
2002) because of their selection. Finally, we compare the 
LF evolution with the predictions of two competitive sce- 
narios of galaxy evolution: the Pure Luminosity Evolution 
(PLE) and the Hierarchical Merging Models (HMM). We 
adopt iJo = 70 km s~^ Mpc~\ = 0.3, Q.a = 0.7. 

2. The K20 spectroscopic sample 

From the total K20 sample of 546 objects with Kg < 20 
we have extracted a sample of 489 galaxies, after exclud- 
ing objects classified as stars and AGN on the basis of 
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redshift 

Fig. 2. K-corrections colors (see text) as a function of red- 
shift: Jrest — Ks (top panel) and K surest — Ks (bottom panel). 
Different curves are derived from different spectral models 
which at low redshift reproduce spectra and colors of local 
E, SO, Sa, Sb, Sc, Im (solid lines from bottom to top at 2 > 1) 
and dusty star- forming galaxies (dashed line; see text). Dotted 
line shows the 2.5log{l + z) term for comparison. 

their spectra. In addition to spectroscopic redshifts, pho- 
tometric redshifts have been derived (Cimatti et al. 2002b) 
with very high accuracy {{zgpe — Zphot) = i'^z) = 0.012 
and CTAz = 0.089{1 + Zspe))- Such an accuracy, made possi- 
ble by the numerous photometric bands and the precision 
of the magnitude measurements, is consistent with previ- 
ous results on the HDF-N data set (Fontana et al 2000) 
and with simulations made by Bolzonella et al. (2000) at 
high redshift (1 < z < 2), where our sample is the largest 
to date. The magnitude-redshift distribution of selected 
galaxies is shown in Fig. ^ Thanks to the depth and com- 
pleteness of our sample, the coverage of the Kg - redshift 
plane is such that the near-IR LF can be derived with 
good accuracy and over a relatively large range of magni- 
tudes at least up to z ~ 1, still sampling with relatively 
good statistics the high luminosity part of the LF up to 
2 ~ 2. Since the sample spans a wide range in luminosity, 
redshift and look-back time, it is therefore well suited to 
study the evolution of the near-IR Luminosity Function 
within the sample itself and in comparison to the local 
population. 

3. The estimate of the rest-frame near-IR 
luminosities 

The small variance of k-corrections in the if^-band (Cowie 
et al. 1994, Mannucci et al. 2001) greatly simplifies 
the analysis with respect to optically-selected samples of 
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redshift 

Fig. 3. The color-redshift {R — Ks vs. z) diagram for galaxies 
in the spectroscopic sample to Ks < 20. Different symbols as 
in Fig. Q Arrows indicate 3a lower limits in R — Ks colors for 
objects not detected in the R band. The curves show different 
spectral models as in Fig. |5| 

galaxies. Moreover, we take advantage of the fact that at 
z ~ 0.8 the observed if^-band corresponds to rest-frame 
J-band. In order to compute the near-IR luminosity func- 
tion, we derive the rest-frame J- and Ks-hand absolute 
magnitudes (Mj and Mk, ) according to the following re- 
lations: 

Mj ^Ks~ 5log{dL{z)/10pc) + {Jrest - K,){z) (1) 
Mk, ^Ks- 5log{dL{z)/10pc) + {Ks,rest - Ks){z), (2) 

where Kg is the total apparent magnitude measured with 
SExtractor (see Cimatti et al. 2002b, Paper III, for more 
details), dL{z) is the luminosity distance at redshift z and 
( Jrest — Ks) and {Ks^rest — Kg) are the "k-correction col- 
ors", i.e. the difference between rest frame and observed 
magnitude which includes also the 2.5Zof/(l -I- z) term (see 
also Lilly et al. 1995). It should be noted that in Eq. |^ the 
last term can be written as a conventional k-correction for 
the Ks-haiid (as in Eq. plus a rest-frame [Z-Ks) color. 
The k-correction colors as a function of redshift are plot- 
ted in Fig. 121 This figure clearly shows the small variation 
of the k-corrections for different spectral types over most 
of the redshift range covered by our data (for a compari- 
son with a similar plot in the optical bands see Lilly et al. 
1995). 

These k-corrections for different spectral types are 
computed using the Bruzual & Chariot (1993) models 
(GISSEL 2000 version), in order to reproduce spectra, and 
k-corrections of local galaxies (E, SO, Sa, Sb, Sc, Im) at 
low redshift (cf. Pozzetti et al. 1996, Mannucci et al. 2001) 
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Fig. 4. The absolute J and Ks magnitudes as a function of 
redshift derived (see text) for all galaxies in the spectroscopic 
sample to Ks < 20. Symbols as in Fig. Local estimates of 
M} and M^, (dotted horizontal lines) from Cole et al. (2001). 

and the color-redshift diagram {R~ Ks vs. z) observed in 
our sample (Fig.lJJ. In order to take into account the dusty 
star-forming population, we introduce an additional spec- 
tral type dominated by star formation and strong dust 
extinction consistently describing the SED of red emis- 
sion line objects and reproducing the sub-class of emis- 
sion line EROs with R ~ Kg > 5 (cf. Paper I). We have 
adopted an exponentially declining star formation history 
{SFR oc exp{-t/T)), with t = 0.3, 1, 2, 4, cx) Gyrs (E, 
SO, Sa, Sb, Sc galaxies, respectively), the Salpeter (1955) 
Initial Mass Function (IMF) and an age of 12.5 Gyrs at 
z = (i.e. a formation redshift zj — 5.7). A model with 
an age of 1 Gyr at all redshifts is adopted for Im galax- 
ies, while a constant SFR and E{B — V) — 0.5 has been 
assumed for the dusty star- forming galaxies (cf. Paper I). 
We have then assigned a spectral model to each galaxy 
according to their R ~ Kg color and spectral type. 

We tested that different assumptions on the k- 
corrections, i.e. by varying model parameters or using 
the fitting of multi-band photometry when available, did 
not affect significantly our results. For example, the k- 
correction estimates which make use of multi-band fitting 
of the SED are in excellent agreement with our estimates 
from R—Kg colors alone, with no systematic shifts in mag- 
nitude and a dispersion of about 0.08 in the J-band and 
0.13 in the K-band up to z 1.5. Varying parameters 
in the synthetic spectral models, i.e. the IMF (Salpeter 
1955, Scalo 1986 and Kennicutt 1983) and the formation 
redshift [zf — 2,3,6), the systematic shifts and disper- 
sions at z < 2 are always less than 0.05 magnitudes. We 
note that magnitudes derived using the k-corrections esti- 



mated from observed spectra of local galaxies (Mannucci 
et al. 2001, Cowie et al. 1994) would be brighter by about 
0.1 magnitudes at z > 0.7, and could be considered up- 
per limits due to the possible evolutionary effect on the 
galaxy spectrum. We are therefore confident that the de- 
rived luminosities are robust with respect to varying this 
set of assumptions. The resulting effects on the Luminosity 
Functions induced by changes in our assumptions are all 
comparable with or smaller than the statistical errors (see 
Sect. El . 

Fig.^shows the absolute magnitudes of all galaxies in 
the spectroscopic sample to Kg < 20 as derived accord- 
ing to our adopted recipes. For the unidentified or unob- 
served galaxies we have adopted the photometric redshift 
Zphot- The photometric threshold of the survey. Kg < 20, is 
clearly visible in the data and translates into different lu- 
minosity limits as a function of redshift, varying from less 
than O.IL* at z < 0.5 to about ^ 0.4L* at z ~ 1, assum- 
ing M} = -23.13 -I- 51og/i7o and M^^ = -24.21 -I- 51og/i7o 
(horizontal lines in Fig. 0)) for the local LF (Cole et al. 
2001). Adopting the mean stellar mass-to- light ratio and 
representative stellar mass {Aigtars/ Lk = 1-32 in solar 
unit and Miters = 1-44 x lO^/ifp^M© for the Salpeter 
IMF) in the local universe (Cole et al. 2001), the lim- 
its in luminosities correspond to 0.1-0.5A^*j^^g or even 
smaller stellar masses since the Aigtars/ Lk increases as 
the stellar population ages (Madau, Pozzetti, & Dickinson 
1998). 

4. The estimate of the near-IR Luminosity 
Function 

After the estimate of the rest-frame near-IR luminosities, 
we compute the J- and Kg-hand Luminosity Functions in 
the following redshift bins: (a) 0.2 < z < 0.65 {zmean — 
0.5), (b) 0.75 < z < 1.3 {z^ean ^ 1-0) and (c) 1.3 < 
z < 1.9 {zjnean — 1-5). The rcdshift bin 0.65<z<0.75 has 
been excluded from this analysis because dominated by 
two groups/clusters (see Fig. 1 and Paper III for details). 
The first two redshift bins include respectively 132 and 170 
galaxies, of which only few objects (2 and 21 respectively) 
do not have a spectroscopic redshift. The highest redshift 
bin includes a smaller number of galaxies (42) with a high 
fraction of photometric redshifts (60%). 

The Luminosity Functions in each redshift bin are es- 
timated using both the 1/Vmax formalism (Schmidt 1968, 
Felten 1976) and the maximum likelihood method STY 
(Sandage, Tammann & Yahil 1979), in order to represent 
our data with a Schechter (1976) parameterization. In the 
^/Vmax analysis, for any given redshift bin (zi, Z2) a max- 
imum volume is assigned to each object. This volume is 
calculated between zi and z„p, the latter being the mini- 
mum between Z2 and Zmax, i-e- the maximum redshift at 
which this galaxy would have satisfied the magnitude limit 
{Kg < 20) of the survey. We use the model tracks in order 
to compute z^ax- Since the STY method determines the 
shape of the LF but not its overall normalization, we have 
normalized the STY LF by matching the galaxy number 
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Table 1. Schechter Luminosity Function parameters ($(M) = 0.4;n(10)$*10-°-''(*^-*^*)("+i)ea;p[-10-°-^(*'^-^^')]) 



Band 


z range 


M - 5 log 


/i70 range 


a 


M* - 5 log hro 


<E>*/i7(f (10-3 Mpc 


J 


0.20 - 0.65 


[-18.7, 


-25.2] 


-1 99+0.22 
-'-•^^-0.20 


OQ QQ + 0.51 
Z,O.Oy_Q gg 


1.99i;* 


J 


0.75 - 1.30 


[-21.8, 


-26.3] 


-0.861^4! 


-23.75t°:« 


3.44tli 


J 


1.30 - 1.90 


[-23.3, 


-26.3] 








Ks 


0.20 - 0.65 


[-19.5, 


-26.0] 


1 9C+0.25 
^•^'-'-0.20 


-24.871^73 


1.781JJ 


Ks 


0.75 - 1.30 


[-22.7, 


-27.2] 


-0.98t»:g 


-24.77l»:S 


2.9lt\l 


Ks 


1.30 - 1.90 


[-24.1, 


-27.1] 









counts observed in each redshift bin. We have checked and 
tested the rehabihty of our results with the independent 
softwares developed by Zucca et al. (1997) and by Poh et 
al. (2001), always finding a good agreement between the 
different methods. 

Because of the relatively small number of objects in 
each redshift bin and only a part of the luminosity func- 
tion can be recovered, the Schechter parameters (a, M* 
and (j)*) derived from our STY analysis are not very well 
constrained by our data. In particular, the uncertainty on 
the faint end slope a increases with redshift because of the 
increase with redshift of the minimum observable luminos- 
ity. A sample reaching ~ 21.0-21.5 would be needed 
to better constrain this parameter at z > 1. The results 
will be discussed in Sectional 

4.1. The treatment of the unidentified objects 

The treatment of the unobserved or spectroscopically 
unidentified sources is not a major problem because of 
the high spectroscopic redshift completeness in our sur- 
vey. Moreover, as discussed in Sect.|2 tested and reliable 
photometric redshifts were derived for most of the faint 
unidentified or unobserved objects, leading to an almost 
negligible number of objects (9 in total) without any red- 
shift information (see Paper III). 

The small uncertainty introduced by the unobserved 
or unidentified objects has been addressed in two different 
ways. 

We first use only the spectroscopic sample and take 
into account the remaining redshift incompleteness (72 
galaxies) by applying weights to each galaxy with spectro- 
scopic redshift both in the Vmax analysis (Avni & Bahcall 
1980) and in the standard STY formulation (Heyl et al. 
1997, cf. also Zucca, Pozzetti & Zamorani 1994). The 
weighting corrections are computed on the basis of the 
fraction of galaxies with spectroscopic redshift in different 
regions of the {R — Kg) — Kg plane. This weighting scheme 
assumes that unidentified objects have the same redshift 
distribution as the spectroscopically identified with simi- 
lar Ks magnitudes and R — Kg colors. While this scheme 
can be reasonably applied over most of the (i? — Kg) — Kg 
plane, the weights are highly uncertain for the optically 
faintest and reddest galaxies, most of which do not have a 



spectroscopic redshift. Since the photometric redshifts of 
these objects are statistically higher than those sampled 
by the measured spectra, this method would not allow to 
correctly estimate the luminosity function in the highest 
redshift bin. 

As a second approach, we therefore base the "best" es- 
timate of the luminosity function on the sample of galaxies 
with spectroscopic plus photometric redshifts. In this case, 
the 9 objects without any redshift information (because of 
the absence of the extended multi-band photometry) have 
an almost negligible effect. To account for them we correct 
the STY normalization (from 1% to 8% in the 3 redshift 
bins) according to the fraction of such objects in each bin. 
The results from the two methods (Fig.|3and|nil are com- 
pletely consistent with each other in the first two redshift 
bins, while some discrepancies are present in the highest 
redshift bin because, as mentioned above, the weighting 
correction function is unable to properly take into account 
the optically faintest and reddest galaxies. Furthermore we 
tested that the uncertainties in the photometric redshifts 
did not affect significantly our results. We extracted ran- 
dom samples from the original photometric redshift cat- 
alog using the measured dispersion. The derived fluctua- 
tions in each magnitude bin of the LFs resulted smaller 
than Poisson errors and in general the simulated LFs are 
consistent with the "best" estimate also in the higher red- 
shift bin with no significant systematic effects. In the fol- 
lowing sections we will discuss the "best" LFs derived us- 
ing spectroscopic plus photometric redshifts. 

5. The cosmic evolution of the near-IR 
Luminosity Function 

Thanks to the high statistical significance and complete- 
ness of our sample it is possible to investigate the evolution 
of the near-IR LF over a wide range of cosmic time. Figs. 
13 and El show the 1/Vmax and the maximum likelihood 
(STY) results for the LFs in the selected redshift bins and 
in the J- and ii's-bands, respectively. 

We find that the luminosity functions are fairly well 
fitted by Schechter functions in the first two redshift bins. 
The magnitude ranges and the best-fit Schechter param- 
eters are summarized in Table ^ with the uncertainties 
derived from the projection of the 68% confidence ellipse. 
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Fig. 5. The rest-frame J-band Luminosity Function in the red- 
shift bins: Zmean — 0.5 (top panel) Zmean — 1-0 (middle panel) 
and Zmean — 1-5 (bottom panel). Points have been derived 
from 1/Vmax analysis (open circles using only spectroscopic z 
and weighted incompleteness corrections, while filled small dots 
using both spectroscopic and photometric z, see Section [4.111 . 
while solid curves are the LF Schechter fits derived from the 
STY maximum likelihood analysis (thin solid lines are the fits 
obtained fixing the a parameter at the local value, see text). 
The dotted curves and vertical dotted lines show the local LF 
in the J-band and M} at « = from Cole et al. (2001). 



In the highest redshift bin, due to the bright and limited 
range in magnitude, our data are well fitted with a sin- 
gle power law (a ~ This does not mean that a 
Schechter function can not be a good representation of the 
luminosity function also in this redshift bin. This is shown 
by the thin solid lines in the lower panels of Figs. and 
which represent fits with a Schecther function obtained 
by fixing the faint end slope a at the local value. However, 
because of the bright limit in absolute magnitude at these 
redshifts, which reduces significantly the range of sampled 
luminosities, the Schechter parameters are very poorly de- 
termined by these data. For this reason the Schechter pa- 
rameters for this redshift bin are not reported in Tabled 
We point out that a direct comparison of the Schechter 
parameters could be quite misleading because the param- 
eters are correlated. 

In order to trace the evolution down to z ~ 0, we 
compared the observed LFs as derived from our sample to 



^-4 
ad 

e-5 



" 9 U 

f^-2 : 



(0 

"-4 
^-5 



o 



-2 
-3 

-4 r 
-5 - 




LLF(K); Cole et aL 2001 
LLF(K); Kochanek et al. 2001 



I I I I 



0.75<z<l 



3 (z =1.0) 

^ mean ^ 




LLF(K): Cole et al. 2001 
LLF(K): Kochanek et al. 2001 



H h 



H 1 



1.3<z<1.9 (z„„=1.5) 

mean / 



LLF(K): Cole et al. 2001 
LLF(K): Kochanek et al. 2001 




-20 



-22 



-24 



-26 



5 log h 



70 



Fig. 6. The same as in Fig.|3 but for the rest-frame Tfs-band. 
The dotted curves and vertical dotted lines show the local LF 
in the Ks-hand and at z = from Cole et al. (2001), while 
dashed curves show the local LF from Kochanek et al. (2001). 

Table 2. Galaxy luminosity evolution at different redshift 



Band 


AM{z = 0.5) 


AM{z = 1.0) 


AM{z = 1.5) 


J 




-0.69«:- 


1 -17+0.22 
^•^ ' -0.24 


Ks 




-0.54«:- 


-1 07+"-^'^ 

l.U( 27 



the Local Luminosity Function (LLF) by Cole et al. (2001) 
in the J and Kg bands (Figs. HI and by Kochanek et 
al. (2001) in the Ks-ha,nd. We then explored two possible 
evolutionary scenarios: luminosity or density evolution. 

Our analysis (both Vmax and STY) shows that the data 
(Figs. Inland are consistent with a mild evolution from 
z = to 2; ~ 1.5 both in J- and Kg-hands. In particular 
the faint-end of the LFs {L < L*) is consistent with the 
local estimates, with no statistically significant evidence 
for a change either in the slope or normalization up to 
z < 1.3 (consistently with Cohen 2002). At higher redshift 
this part of the luminosity function is not well sampled by 
our data, but also in this bin a Schechter function with 
local faint-end slope is consistent with the data (Figs. El 
andEIl. Viceversa, it is interesting to note that the density 
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Fig. 7. Rest-frame Ks-hand Luminosity Function in two red- 
shift bins: Zmean — 0.5 (top panel) Zmean — 1-0 (bottom panel). 
Points and lines derive from 1/Vmax and STY analysis for early 
(filled symbols and dotted lines) and late (open symbols and 
dashed lines) galaxies respectively. The solid curves are the 
STY LF for the total sample in the Ks-hand as in Fig. [5] 



of luminous galaxies (e.g. Mk^ — 5 log/170 < —25.5) is 
significantly higher than locally. 

Because of the statistical correlation between the 
Schechter parameters {a, M* and $*), in order to estimate 
the luminosity and/or density evolution within our sample 
and in comparison with the local LFs, we have fixed the a 
parameters of our LFs to the same values observed locally 
(Cole et al. 2001), a,/ = -0.93 and a^, = -0.96 (thin 
solid lines in Figs. [3 EJ. We are allowed to do this be- 
cause all our best fit a values are consistent with the local 
ones (see Table^l. Assuming local M} = — 23.13-|-51og/i7o 
and M^^ = -24.21 + 51ogft,7o from Cole et al. (2001) 
we find a luminosity evolution at 2 ~ 1 of the order of 
AMj ~ -0.69 ± 0.12 and AM^ ^ -0.54 ± 0.12 with 
the normalizations (0*) in the first two redshift bins con- 
sistent with the local values, within from 3 to 20%. At 
z > 1.3 there is an indication of an even higher luminos- 
ity evolution (AMj ~ -1.17, AM/^ ~ -1.07), while the 
(p* values decrease by a factor 4-5 with respect to the lo- 
cal value. This last result should be taken with caution 
because of the lower statistics in the highest redshift bin 
and the smaller range in luminosity covered. 

The luminosity evolution estimated from STY analysis 
are summarized in Table |21 with the formal statistical la 
confidence limits. The uncertainties introduced by the k- 
correction, discussed in Section|31 are always less than the 
statistical errors, and less than 0.03 and 0.10 magnitudes 
at 2; < 1.3 and z ~ 1.5, respectively. 



We recall here that Cimatti et al. (2002b) discussed 
the photometric selection effects present in the K20 sam- 
ple and showed that, on average, the total flux of spirals 
and ellipticals with L < L* are underestimated by 0.1 
and 0.25 magnitudes respectively, while at higher L the 
flux lost for elliptical galaxies could be even higher. For 
this reason we are confident that our estimate of a posi- 
tive luminosity evolution is quite conservative. The above 
estimates of the luminosity evolution could increase on av- 
erage by about —0.2 magnitudes if we take into account 
this underestimate of the total flux. 

6. The Luminosity Function by spectral or color 
type 

In order to investigate the role of different galaxy pop- 
ulations, we divide the K20 sample in two subsamples 
on the basis of the spectroscopic classification, i.e. early 
type galaxies without strong emission lines and late type 
galaxies with emission lines, and study the near-IR LF 
for the different populations in the two redshift bins at 
Zmean — 0.5 and 1. The high fraction of photometric red- 
shifts without any spectroscopic information does not al- 
low us to extend such studies to the highest redshift bin. 
We take into account the spectral incompleteness using 
two weighting functions (see Sect. I4.1|l derived for each 
spectral type and correcting the overall normalization in 
the STY method using the color information of unclassi- 
fied objects. 

The results are shown in Fig. [3 Our analysis indi- 
cates that the faint-end slope of the LFs at these red- 
shifts shows the same general dependence on galaxy spec- 
tral type as that found in numerous studies of the local 
LF. Specifically, the faint-end slope is much steeper for 
late/emission line galaxies than for early type galaxies. 
A similar result has been found by Cohen (2002) in the 
HDF-north. More important, we find, for the first time, 
that at Zmean — 0.5, carly type galaxies clearly dominate 
the bright-end of the Luminosity Function, in agreement 
with local observations (Kochanek et al. 2001). A similar 
result is also visible at z^ean — 1-0. Some of the pre- 
vious surveys, selected in the optical bands with follow- 
up if-band observations (e.g. Cohen 2002), present a bias 
against identifying red early type galaxies at high redshift 
because of substantial incompleteness of EROs in the sam- 
ples. 

The overall density of early type galaxies at 2 ~ 0.5 is 
consistent with the local estimate (Kochanek et al. 2001), 
while at z ~ 1, their number density shows at most a small 
decrease (< 30%), consistent with recent results from Im 
et al. (2002). This decrement could be due both to the 
existence of a population of galaxies morphologically clas- 
sified as early-type but with blue colors due to some recent 
episodes of star formation (Menanteau, Abraham & Ellis 
2001, Im et al. 2001) and to our spectral incompleteness, 
more severe for red/early type galaxies. 

Similar results are obtained if we divide the galaxies 
in two samples according to their R — Kg colors (we adopt 
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Fig. 8. — ATs colors vs. rest-frame absolute ifs magnitudes 
for K20 galaxy sample in 3 different redshift bins. Circles rep- 
resent galaxies identified spectroscopically, while squares are 
unidentified or unobserved galaxies plotted a,t z = Zphot- In 
each redshift bin empty and filled symbols refer to galaxies 
with z < Zmean and z > Zmean, respectively. 



the color of Sa galaxies to divide the sample in red and 
blue). 

We conclude that red and early type galaxies dominate 
the bright-end of the Luminosity Function already at early 
epoch and that their number density shows at most a small 
decrease (< 30%) up to z 1. 

Fig- IHl shows the optical/near-IR colors versus near-IR 
luminosities in the three redshift bins. While the magni- 
tude limit of the survey, i^s < 20, corresponds to different 
low luminosity limits in the different redshift bins, there 
appears to be a correlation of the optical/near-IR colors 
with near-IR luminosities. In first approximation this cor- 
relation can be interpreted as a correlation of age and/or 
specific star formation rate with stellar mass, with the 
most "massive" galaxies being "old" and the "low-mass" 
galaxies being instead dominated by young stellar pop- 
ulations. This is now well established for local galaxies 
(Gavazzi et al. 1996, Boselli et al. 2001, Kauffmann et 
al. 2002) and we extend it to z ^ 1.5. Fig. |S1 (see empty 
and filled circles representing galaxies with z < Zmean and 
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Fig. 9. Comoving luminosity density in J— (top panel) and 
7('s-bands (bottom panel). The filled circles are the value de- 
rived directly from observations (see text), while open circles 
are the "LF-corrected" estimates. The solid lines show the best- 
fit power laws. The open squares are from Cole et al. (2001). 



duced by a redshift effect within the bins. In Sect. |Slwe 
compare the observed correlation with HMM predictions. 

7. The near-IR luminosity density evolution 

Tracing the cosmic emission history of the galaxies at dif- 
ferent wavelengths offers the prospect of an empirical de- 
termination of the global evolution of the galaxy popu- 
lation. Indeed, it is independent of the details of galaxy 
evolution and depends mainly on the star formation his- 
tory of the universe (Lilly et al. 1996, Madau, Pozzetti 
& Dickinson 1998). Various attempts to reconstruct the 
cosmic evolution of the comoving luminosity density have 
been made previously mainly in the UV and optical bands 
(Lilly et al. 1996, Cowie et al. 1999). Our survey offers the 
possibility to investigate it in the near-IR using a LF ex- 
tended over a wide range in luminosity. 
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We estimate the comoving luminosity density directly 
from observed sources using the Vmax formahsm (see for 
details Lilly et al. 1996). The major source of uncertainty 
is the contribution from galaxies fainter than the magni- 
tude Hmit Kg < 20. We have therefore estimated also the 
"LF-corrected" luminosity density using the best-fit LF 
from STY method and extending it to fainter luminos- 
ity (we adopt 0.02L*(z) as our lower limit in luminosity 
for this extrapolation) . A formal uncertainty in this proce- 
dure was estimated by considering the range of acceptable 
Schechter parameters values (see Table P). 

Fig. 1^1 shows luminosity densities derived in J- and 
i^s-band as a function of redshift. Up to z < 1.3 the un- 
certainties due to faint galaxies and LF parameters are 
small. At higher redshift the uncertainties in the LF pa- 
rameters do not allow to constrain the luminosity density 
values. We could therefore only estimate a lower limit to 
the luminosity density using the observed value. 

We find a slow evolution with redshift of the observed 
near-IR luminosity density, consistent with the results by 
Cohen (2002) in the HDF-north, even if she had to apply a 
large correction to her data because of the disappearance 
of absorption galaxies at high redshift in her sample as 
they became EROs. Using the observed local luminosity 
densities derived from Cole et al. (2001), the luminosity 
density evolution up to z < 1.3 is well represented by a 
power law, pxiz) = px{z = 0)(1 + z)'^(^). We find (3{J) ~ 
0.70 and ^{Ks) ~ 0.37 (Fig.|3. 

The near-IR luminosity density evolution is much 
slower than that found in the UV and optical bands 
(/3 = 3.9-2.7 from 0.28 to 0.44 /im by Lilly et al. 1996 
and (3 = 1.5 at 0.15, 0.28 ^m by Cowie et al. 1999, 
for = !)• Indeed, while the optical luminosity den- 
sity evolution is mainly related to the star formation his- 
tory, the evolution of the near-IR luminosity density is 
more closely related to the stellar-mass density. The lo- 
cal stellar mass density, derived by Cole et al. (2001) is 
fl^tars = (3.7 ± 0.6) X 10^^/ifo^. If we adopt the stellar 
mass-to-light ratio for the galaxy spectral models which 
best match colors and spectral types for each galaxy, the 
mean Adgtars/ Lk in our sample becomes 0.63 and 0.54 
at z = 0.5 and 1, with small variations due to the model 
parameters adopted (see Sect. O. With these values of 
Mstars/Lx we derive flstars = (2.0 ± 0.1) x lO^^/ifo^ 
and ^stars = (2.1 ± 0.3) x lO'^hj^ at z = 0.5 and 1 
respectively. This analysis suggests that the evolution of 
the stellar mass density is relatively slow with redshift, 
with a decrease of about a factor 1.8 ± 0.4 from z = to 
z ~ 1, consistently with recent results from Dickinson et 
al. (2002a). 

At z > 1 we can only give a lower limit to the near- 
IR luminosity density. As discussed by Madau, Pozzetti, 
Dickinson (1998) (cf. also Pozzetti & Madau 2001), a ro- 
bust determination of the ncar-IR luminosity density at 
z > 1 could be fundamental to disentangle between dif- 
ferent cosmic histories of star formation (SFH). In fact, 
SFHs peaked at intermediate redshift (1 < z < 2) predict 
a decrease in the near-IR luminosity density at z > 1-2, 



while if 50% of the stars formed at z > 2 (similar to the 
PLE adopted here), the corresponding near-IR luminosity 
density flattens at z > 1 without a strong decrease (cf. 
Fig. 6 in Pozzetti & Madau 2001). A detailed comparison 
of luminosity and stellar-mass density with model predic- 
tions and empirical star formation histories will be further 
investigated in a forthcoming paper. 

8. The comparison with model predictions 

We have compared our LF with the predictions from two 
competitive scenarios of galaxy evolution: the Hierarchical 
Merging Models (HMMs) and the Pure Luminosity 
Evolution (PLE) (Fig. Unj. For the HMMs, we used 
the predictions by Cole et al. (2000, COO, kindly pro- 
vided us by C. Baugh), Menci et al. (2002, M02), and 
GIF (Kauffmann et al. 1999). For the GIF simulations, 
which combine the large high-resolution N-body simula- 
tions with semi-analytical models, we derive the LF at a 
magnitude above which it can be considered complete in 
luminosity given the mass limit of the simulated galaxy 
sample (2 x 10^° H^^Mq). For the PLE parameterization, 
we adopted the predictions of Pozzetti et al. (1996, 1998, 
PPLE) (for details see also Paper IV). 

For a more reliable comparison with models, we cor- 
rected the observed magnitudes by —0.2 mag, in order to 
take into account the average loss of total flux in the K20 
sample (see Papers HI- IV). 

As discussed in Paper IV at low-z (z — 0.5) the HMMs 
by Cole et al (2000) and by Menci et al. (2002) predict 
many more galaxies than observed (cf. Fig. 3 in Paper 
IV) . From the present study (upper panel of Fig. I10|l we 
can therefore conclude that such excess is due to low lu- 
minosity, "low-mass" objects, which dominate the steep 
faint-end in the HMM LF. At z = 1, the HMMs predict 
even steeper LFs because low-mass galaxies become more 
numerous in the hierarchical scenario at higher redshift. 
This is in contrast with the LF derived from our data at 
z = 1, which does not show a similar behaviour. These 
discrepancies are consistent with previous results both in 
the local universe (cf. Fig. 1 in Baugh et al. 2002) and at 
high redshift (z < 1 and z 3) for the B and UV rest- 
frame LF, respectively (Poli et al. 2001, Somerville et al. 
2001). 

Also at the bright end of the LF the HMMs appear to 
be in disagreement with our data at z ~ 1, where they un- 
derpredict the density of bright {Mk^ < — 25.5 + 51og/i7o) 
galaxies with respect to our data. For example, our data 
at Mk^ = —25.6 + BloghjQ and z ~ 1 are about a fac- 
tor 2.6 higher than the HMM predictions and even higher 
at brighter magnitudes if we consider the STY maximum 
likelihood fit. In the highest redshift bin the excess over the 
predictions becomes significant at AIk^ < —26-1-5 loghjQ. 
In comparison, the GIF simulated catalogue at z — 1 un- 
derpredicts the K20 LFs at all magnitudes. As discussed 
by Kauffmann et al. (1999), the GIF model produces a 
factor 2-3 too few galaxies at magnitudes around L* al- 
ready at z = 0, while it produces an excess of bright local 
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Fig. 10. The rest-frame Ka-hand Luminosity Function in the redshift bins; Zmean — 0.5 (top panels) Zmean — 1-0 (middle 
panels) and Zmean — 1-5 (bottom panels) compared to PLE (right panels) and HMM (left panels) models. We have used 
the models closer to Zmean or interpolated between different redshifts. Data points derive from 1/Vmax analysis, while dotted 
curves are the LF Schechter fit derived from STY analysis. We have corrected the data on average by —0.2 magnitude because 
photometric selection effects (see text). 



galaxies. This problem -with respect to the local LF affects 
also the Cole et al. (2000) model (see Fig. 4 in Baugh et al. 
2002). Therefore, as discussed in Sect.[5l -while the density 
of luminous objects {Mk^ —5 loghjo < —25.5) is quite con- 
stant or mildly increasing ■with z within our sample, and 
higher than the local density, it is rapidly decreasing in 
the Cole et al. (2000) and GIF models at z > 0.5 (Fig. [TT|l 
in clear conflict -with our data. The problem of a negative 
evolution with redshift of luminous objects does not affect 
instead the Menci HMM model, even if also this model sig- 
nificantly and systematically underpredicts their number 
density at all redshifts. 

In addition to the comparison of the data and predic- 
tions for the LF, it is also important to verify whether the 
bright L > L* galaxies in the K20 survey have colors con- 
sistent with the HMM predictions. Fig. 1 121 shows the com- 
parison between the z — 1.05 GIF simulated catalogue and 
data at 0.75 < 2 < 1.3 {zmean — !)• While the deficiency 
of simulated blue low-luminosity {R — Kg < 4.2, Mk^ — 5 
log /i7o > —23.5) galaxies is due mainly to the mass limit 
of the GIF catalogue, a serious discrepancy emerges in the 
two distributions in the magnitude range where the GIF 
catalogue is expected to be complete. The GIF simulated 
catalogue shows an excess of blue high-luminosity galaxies 
{R-Ks < 4.4, Mk, - 5 log hjo < -24.5) and a deficiency 
of red luminous galaxies (see also the color distributions of 
high-lmninosity galaxies, Mk^ — 5 log hjo < —24.5, in the 
right panel of Fig. 1121 normalized to the same comoving 
volume). This result is in agreement with the finding of a 



deficiency of EROs in the HMMs (Daddi et al. 2000, Firth 
et al. 2002, Smith et al. 2002), but, since it is based not 
only on colors but also on the luminosity of the galaxies, 
it points out in particular a clear deficiency of red lumi- 
nous galaxies in the HMM predictions. This plot clearly 
shows that a comparison of LF and colors is much more 
powerful than that of the colors or LF alone and that 
the GIF Ks-ha,nd LF at z = 1 is obtained with a radi- 
cally different mix of galaxies compared to the observed 
one. While the bright end of the HMM LF at z = 1 is 
dominated by actively starforming galaxies, in the K20 
sample it is dominated by the passively evolving galax- 
ies. Given that blue galaxies suggest younger stellar pop- 
ulation and lower Adstars/ Lk ratios, we expect that the 
failure of this model in reproducing our LF (see Fig. I10|l 
should be even stronger in term of density of "massive" 
galaxies at z 1 (cf. also the luminosity and stellar-mass 
functions in Baugh et al. 2002). 

The right panels of Fig.llOlshow the comparison of our 
data with the PLE predictions, computed with two dif- 
ferent IMF (Scalo and Salpeter). The overall agreement 
is much better than with the HMM predictions even if 
some discrepancies, but at much lower level of significance, 
are present. In particular, in the highest redshift bin the 
parameterization with the Salpeter IMF overpredicts the 
total number of galaxies (see also Paper IV), while the 
Scalo IMF (dominated by low/intermediate mass stars, 
which induce a lower luminosity evolution compared to 
Salpeter at z > 1) fits well the density of galaxies at 
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Fig. 11. Comoving number density of luminous galaxies 
{Mks < —25.5 -I- 5 log/i7o) in the J^s-bands. The open cir- 
cles are the values derived directly from observations, while 
the filled square is from Cole et al. (2001). The different lines 
show models as in figure 1101 
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Fig. 12. Left panel: R — Ks colors vs. rest-frame absolute Ks 
magnitudes for z = 1.05 GIF simulated catalogue (small dots) 
and data (circles) at 0.75 < z < 1.3 {zmean — 1) (empty and 
filled circles refer to 2 < 1 and z > 1 respectively). Vertical 
dashed line represents approximately the completeness mag- 
nitude limit of GIF catalogue corresponding to its mass limit 
(see text). Right panel: Color distribution of high-luminosity 
galaxies {Mk^ — 5 log /170 < —24.5) observed (dotted line) and 
simulated (continue line), normalized to the same comoving 
volume. 



Mk^ > —26.0 -I- 5 log/170, but falls below the data for the 
brighter galaxies. In Paper IV we found that in the con- 
text of PLE models also the observed redshift distribution 
of K20 galaxies is better fitted assuming an IMF dom- 
inated by low/intermediate mass stars (e.g. Scalo IMF) 
than a Salpeter IMF (see also Pozzetti et al. 1996 and 
Broadhurst & Bouwens 2000). PLE models with a flatter 
IMF (i.e. Salpeter-like) can be made consistent with the 
data (redshift distribution and LF) only if the rapid evo- 
lution induced by massive stars is sufficiently suppressed 
by dust attenuation (Totani et al. 2001, see Paper IV). 
Indeed, such a dusty phase is predicted by many models 
for the formation of high redshift spheroids and could be 
associated to high redshift SCUBA sources. 

9. Summary and discussion 

The cosmic evolution of the iiT^-band selected field galaxy 
population has been studied over the redshift interval 
0.2 < z < 1.9 using about 500 galaxies from the K20 
spectroscopic survey. The sample spans a wide range in 
redshift and look-back time and allows to study the evo- 
lution of the near-IR Luminosity Function both within the 
sample and in comparison to the local population. 

We take advantage of the near-IR selection (in partic- 
ular in the iiT- band), in which the k-corrections are rela- 
tively invariant to galaxy type and relatively small also at 
high redshift (Cowie et al. 1994), and the dust extinction 
effects are less severe than in optical samples. 

We derived the near-IR luminosity function in the rest- 
frame J and i^s-band in three redshift bins (zmean — 
0.5, 1.0, 1.5). The detailed analysis of the observed LF at 
different redshifts and the comparison with the local LF 
and with the predictions of galaxy formation models pro- 
vided the following results: 

1- A mild evolution is observed both in the J and Kg 
Luminosity Functions to z ~ 1.5, in agreement with previ- 
ous indications by Cowie et al. (1996), Cohen (2002) and 
Feulner et al. (2003). There is no evidence of a steepen- 
ing of the faint-end LFs up to 2; ~ 1.3. In particular, the 
faint-end (L < L*) is consistent with the local estimates 
up to z < 1.3, while at the bright end the data show a 
luminosity evolution of about AMj ~ —0.69 ±0.12 and 
AMk ^ -0.54 ± 0.12 at z ~ 1. The density of luminous 
galaxies {Mk^ — 5 log/170 < —25.5) is relatively constant 
or mildly increasing with redshift within our sample and 
higher than locally at all redshifts. 

2- Pure density evolution cannot reproduce the ob- 
served LF at z < 1.3. 

3- Red and early-type galaxies clearly dominate the 
bright-end of the Luminosity Function at z ~ 0.5 and a 
similar trend is also visible at z ~ 1, showing that such 
systems were already in place and fully assembled at that 
cosmic epoch. Their number density shows at most a small 
(< 30%) decrease up to z ~ 1. 

4- The evolution of the rest-frame near-IR comoving 
luminosity densities up to z ~ 1 can be described by power 
laws, px{z) = px{z = 0)(1 + zf^^\ with (3{J) ~ 0.70 and 
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(3{Ks) — 0.37. Such an evolution is much slower than those 
observed in the UV and optical bands. 

5- The hierarchical merging models overpredict the LF 
at low luminosity at all redshifts, whereas they under- 
predict the density of high-luminosity galaxies at z > 1. 
HMMs by Kauffman et al. (1999) and Cole et al. (2000) 
overpredict high luminosity galaxies at z = 0, and predict 
a negative density evolution of the bright end of the LF 
at z > 0.5 which is not observed. 

6- The PLE predictions are in rather good agreement 
with the mild luminosity evolution observed up to z ~ 1.5. 

7- There appears to be a clear correlation of the 
optical/near-IR colors, i.e. in first approximation the spe- 
cific star formation rate and ages, with near-IR luminosi- 
ties, i.e. the stellar mass. This correlation suggests that the 
most "massive" galaxies are "old", while the "low-mass" 
galaxies are instead dominated by young stellar popula- 
tions. The GIF model shows instead a clear deficiency of 
red luminous galaxies at z 1 compared to the observa- 
tions. 

From the analysis of the LF and of the comoving 
near-IR luminosity density, we can derive some indica- 
tions on the evolution of the Galaxy Stellar-Mass Function 
(GSMF). A detailed analysis will be presented by Fontana 
et al. (2003). The present analysis suggests that the evo- 
lution of the GSMF and of the stellar mass density is slow 
with redshift up to z ~ 1.5, in contrast with the rapid evo- 
lution of the Galaxy Stellar Mass Function expected in the 
hierarchical models at z < 2 (cf. Baugh et al. 2002). In our 
sample, using galaxy spectral models which match colors 
and spectral types, we found that flstars decreases by a 
factor about 1.8 ± 0.4 from z = to z ~ 1, mainly due to 
the decrease in the mass-to-light ratio, which is however 
slower than in HMMs models. In at least one of these mod- 
els (GIF) the small adopted M/L ratio corresponds to a 
color distribution of galaxies at z ~ 1 clearly inconsistent 
with data (see Fig. EJ- 

Moreover, the fact that the bright-end of the 
Luminosity Function is dominated by red/early type 
galaxies and that the mean spectra of early type EROs 
at z ~ 1 are consistent with an old stellar population (see 
Paper I), suggest that old and massive elliptical galaxies 
were already in place at z ~ 1 (see also Paper I) and, there- 
fore, they should have formed their stars and assembled 
their mass at higher redshift. These results are in contrast 
with the current renditions of hierarchical models in which 
the bulk of massive elliptical galaxies forms through merg- 
ing at low redshift (z < 1-2) and suggest that most of the 
merging which form elliptical galaxies in the hierarchical 
models should occur at higher redshift (say, z > 2-3), fol- 
lowed by a pure "luminosity-like" evolution. On the other 
hand, HMMs overpredict the density of low-luminosity, 
low mass, galaxies at z < 1, as previously noted in the 
local universe (cf. Fig. 1 in Baugh et al. 2002) and in the 
B and UV rest frame bands at z > 1 (Poll et al. 2001, 
Somerville et al. 2001, Poh et al. 2003). 

Our results show the importance of sampling the faint- 
end of the LF in even deeper near-IR selected samples. 
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Fig. 13. Predicted SIRTF-flRAC number densities for galax- 
ies with z < 2 derived from K20 sample (see text) as a 
function of AB magnitudes at 3 different wavelengths (A — 
3.6, 4.5, 5.8/im). Histograms refer to Ks < 20 galaxies in our 
sample in the indicated redshift ranges, while lines have been 
derived from K20 LFs extended to fainter luminosities (using 
the fiat slope in the highest z range). Vertical dashed lines refer 
to GOODS flux limits at the different wavelengths. 



In addition, to derive the LF at z > 1 in a wider 
but shallow sample (say K < 18.5) will be essential to 
trace the evolution of the most massive galaxies with 
higher statistical significance. In the near future direct 
measurements of rest-frame K-hand magnitude for high 
redshift galaxies will be possible with space-based in- 
frared surveys, in particular the SIRTF (Space Infrared 
Telescope Facility) Legacy Science project GOODS (Great 
Observatories Origins Deep Survey, P.I. M. Dickinson). 
Our data allow to directly estimate the number densities 
of z < 2 galaxies at the different SIRTF-GOODS wave- 
lengths, A — 3.6, 4.5, 5.8/im, which sample the rest-frame 
i^s-band at z ~ 0.7,1.1,1.7 respectively, similar to the 
Zmean adopted for the K20 LFs. We have converted the ab- 
solute magnitudes, derived here for the K20 sample, 
to SIRTF-f IRAC fluxes using the galaxy spectral models 
described in Sect. |31 for normal galaxies and the spectral 
model for M82 (Silva et al. 1998) for the dusty star form- 
ing galaxies. The expected number densities for objects 
with z < 2 are shown in Fig.^| At the expected depth of 
GOODS (see Dickinson et al. 2002b), SIRTF will detect 
about 5000 galaxies at z < 2 within 330 arcmin^. GOODS 
SIRTF Legacy Program, sampling the rest frame near- 
IR luminosities directly, will allow to derive the i^-band 
Luminosity Function along with the stellar mass distri- 
bution well below the K20 limit (see Fig. I13|l and with 
higher statistical significance. With these data the stellar 
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mass assembly history of galaxies will be measured over a 
wide range of redshifts and cosmic time. 

We make publicy available the 1/Vmax estimates at 
|http: / / www.arcetri.astro.it~k20 / releases/index.html, 
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